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variant Hamiltonians and study the time average expectation values and entanglement entropies in 
finite chains. We show that in noninteracting models the time average of spin correlation functions 
is asymptotically equal to the infinite time limit in the infinite chain, which is known to be described 
by a generalized Gibbs ensemble. The equivalence breaks down considering nonlocal operators, and 
we establish that this can be traced back to the existence of conservation laws common to the 
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correction for time average correlation functions and entanglement entropies. We find that large 
corrections are generally associated to observables with slow relaxation dynamics. 
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I. INTRODUCTION 

The experimental realization of systems that evolve 
unitarily for long time without dissipation 1 " 5 has spurred 
the theoretical research on nonequilibrium dynamics in 
isolated many body quantum systems. Trapped ultra- 
cold atom gases have become the paradigm of systems 
that are weakly coupled with the environment and can 
be designed to have highly tunable Hamiltonian param- 
eters!^. At the same time, the evolution after a sudden 
quench of a Hamiltonian parameter has come to be one of 
the most studied protocols^!! primarily but not exclu- 
sively because of the stationary properties that emerge 
at late times after the quench. Although the system is 
in a pure state, generally correlation functions as well 
as reduced density matrices (RMDs) relax to stationary 
values^ which can be described in terms of an effective 
statistical ensemble. In addition, the stationary state 
turns out to be strongly influenced by dimensionality and 
local conservations laws^l. 

In this work we consider the time evolution of the 
ground state of a translational invariant spin Hamilto- 
nian after a sudden quench of a global parameter ho — > h 

|*(t)> = e- iH W |* > H[ho] |* ) = E GS |* > • 

The late time physics can be investigated in alternative 
ways. In infinite chains we can study the behavior of 
RDMs at long time after the quench, and define the sta- 
tionary state poo as follows 



lim 



lim 

oo t— >oo L 



lim Tr A [p( L \t)] 



(1.1) 



where p (i) (*) = is the density matrix in a 

chain of L spins, A is the complement of the subsystem 
A, and we assumed the limits to exist (but there are 
exceptions^, e.g. if there are localized excitations). 

Because of quantum recurrence^! p^ cannot be defined 
in finite systems. However generally it turns out that 
at almost every time expectation values of observables 
are close to their time averages^!. Thus, in a sense, the 



finite system relaxes to the ensemble p that describes 
time average expectation values 



p = lim p( L ' 

L— f oo 



Tr[Op (i) ] = lim - / dtTr[Op (i) (i)] VO . 



(1.2) 



It is commonly called diagonal ensembl d 30 l 34 ^ i, since in 
a basis that diagonalizes the Hamiltonian only the di- 
agonal elements of the density matrix survive the time 
average. 

The thermodynamic limit and the infinite time limit 
in eqs ( |1.1[ )( |1.2| come in different order, so a preliminary 
question is whether the two ensembles are identical 

P = Poo- (1.3) 

In general the answer is no. The breakdown of relations 
like (1.3 1 in the nonequilibrium evolution of a generic 
pure state is a well-known factP^ that can be traced back 
to the issue of ergodicity in quantum mechanics^. 

In this paper we focus on one aspect that rules out 
Eq. (1.3). We indeed show that it is not unusual that 



the Hamiltonians before and after the quench share lo- 
cal conservation laws, and this is sufficient to contradict 



Eq. (1.3) 



A second question is whether the two ensembles have 
the same local properties, i.e. 



lim Tr A [p] = p a 



(1.4) 



In other words we wonder whether the expectation value 
of local observables in finite chains are equal, on time av- 
erage, to the va lue t hey approach in the limit 1 <C t <C L. 
We regard Eq. ( 1.4 ) as a nontrivial statement, though the 



physical intuition suggests (1.4) to be true. For a generic 



system, Eq. ( 1.4) together with the cigenstate thermal- 
ization hypothesi s ! 39 * 40 ! underlie the idea of thermal re- 
laxation 41 , which results in an equilibrium thermal en- 
semble poo- It is well known that local conservation laws 
after the quench invalidate this simple description, and 
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it was proposed in RefP^ that the constraints of conser- 
vation laws Ij can be taken into account by considering 
the generalized Gibbs ensemble (GGE) 



GGE 



expt-^A^-), with [!<,!,] =0, (1.5) 
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where I\ = H \h\. T his conjecture has been verified in 
several systemg^EUMl anc [ [ s now expected to describe 
the late time behavior after a quantum quench in inte- 
grable models, in which there are infinite local conser- 
vation laws. Although the GGE 'solves' the problem of 
ergodicity at infinite time after the quench for local de- 
grees of freedom, the fact that poo can be different from 
p shows t hat there are other effects. 

If Eq. (1.4 1 is satisfied then the question moves to the 
finite-size corrections of p^ in Eq. (1.2) and to the ap- 



proach to the stationary state in the infinite chain 
(see also Ref d 35 * 47 ^^). In order to investigate these ques- 
tions we consider quenches in models with a free-fermion 
representation, in which the late time behavior of correla- 
tion functions and entanglement entropies is known, and 
we compute the leading finite-size correction for time av- 
erage correlation functions and entropies. We recognize 
that large finite-size corrections can be ascribed to the 
existence of (local) conservation laws 'unaffected' by the 
quench. Moreover, we find that they are generally asso- 
ciated to observables with slow relaxation dynamics. 



Quenches in noninteracting spin chains: 
the price of solvability 

Noninteracting spin chains are the standard testing 
ground in noncquilibrium problems. If at any time the 
Hamiltonian can be mapped into a quadratic fermionic 
operator by a time independent transformation, the sys- 
tem can be studied by means of standard free fermion 
techniques. This allows numerical analyses in polyno- 
mial computational time, as well as analytic investiga- 
tions. We point up that most of the simplifications rely 
on the time independence of the fermion mapping. 

We focus on periodic spin chains (including chains that 
can be interpreted as translational invariant spin ladders) 
and call ag the Majorana fermions in terms of which the 
Hamiltonian is quadratic. If the mapping is given by 
a Jordan- Wigner transformation the Majorana fermions 
can be defined as follows 



an = Y[ 



a 2 *-!=n^ (1-6) 



and we have {cii,aj} = 26ij. Short range and transla- 
tional invariance result in a Hamiltonian of the form 



(1.7) 



' 71=1 



with % a block circulant matrix (apart from boundary 



matrices (i, j — 1, . . . , L/d, in contrast to the indices in 
Eq. (1.7)). The dimension d of the blocks depends on 
the spin interaction. On the other hand the circulant 
structure is a consequence of translational invariance and 
hence is independent of the Hamiltonian parameters. We 
note that relaxing the condition of short range is not 
generally sufficient to spoil the structure. 

It is straightforward to show that in the infinite chain 
there are infinite conservations laws independent of the 
system details. Indeed, if C is a generic L/d circulant 
matrix we have 

[- a Me n a n , | at[C ® I d ] 

t,n t,n 

where in the last step we used that circulant matrices 
commute between each others. This shows that in short- 
range noninteracting Hamiltonians translational invari- 
ance produces infinite conservation laws independent of 
the Hamiltonian parameters. Moreover, their range is es- 
sentially given by the size of the block of diagonals around 
the main diagonal that include all nonzero elements of C, 
so they can be defined to have finite range. 

The dynamics can be solved because it reduces in fact 
to a d-dimensional problem, and hence implicitly because 
of the infinite local integrals of motion independent of 
the system details. Although quenches in noninteracting 
models have been widely investigated, not much empha- 
sis has been placed on this point and, to our knowledge, 
no observed behavior has been ascribed to this very as- 
pect. We are going to show that its signature can be 
found in the finite-size corrections and in the relaxation 
dynamics of observables. 

The manuscript is organized as follows. In Section [IT] 
we introduce the model and reportP^ the explicit ex- 
pressions of Eqs (1.1)(1.2). In Section III we show that 



the diagonal ensemble in the transverse-field Ising chain 
(TFIC) describes noninteracting fermions and we con- 
struct the noninteracting representation; we also discuss 
the effects of common conservation laws. Section IIVI is 



quite technical: we prove Eq. (1.4) in the TFIC and de- 



terms), that is to say Hi 



Ha- 



(i—j) mod (L/d) 



are d x d 



velop a formalism to compute the leading finite-size cor- 
rection of correlation functions. In Section [V] we study 
the finite-size corrections for the two-point functions of 
the most important operators of the model. We also com- 
pute the finite-size correction for the Renyi entropy S2 
and discuss the correction for the von Neumann entropy. 



II. THE MODEL 

In order to draw a comparison between relaxation 
to the stationary state in the thermodynamic limit 
and finite-size corrections in finite chains we focus on 
quenches in the transverse field Ising chain, where re- 
cently the time dependence of correlation functions has 



3 



been computed exactly^ 28 ' 51 1 53 1 54 ! We note that our dis- 
cussion generalizes straightforwardly to other spin chains 
with free fermion spectra, like the quantum XY model. 
The Hamiltonian is given by 



(2.1) 



At zero temperature the model exhibits ferromagnetic 
(h < 1) and paramagnetic (h > 1) phases separated 
by a quantum critical point h c = 1. It is the simplest 
paradigm of quantum critical behavior and quantum 
phase transitions, and in the last years has also become 
a crucial paradigm of quench dynamics. From a tech- 
nical point of view, the Hamiltonian (2.1) with periodic 



boundary conditions is mapped into two separate nonin- 
teracting fermionic sectors by the Jordan- Wigner trans- 
formation (1.6) and is finally diagonaliz ed by a Bogoli- 
ubov transformation in momentum space^^ (d = 2). 

In the finite chain the state before the quench, namely 
the ground state of the model with a given magnetic field 
h , is the squeezed coherent state^ 



I* 



o ) oc exp 1 1 



E 

0<pSNS 



K(p)a 



t t 
a. 



10) 



(2.2) 



where |0) is the vacuum of the a fermions that diagonal- 
ize the antiperiodic (Neveu-Schwarz) sector of the final 
Hamiltonian as follows 



Hns = 



J2 Sh 

pGNS 



0) (<4< 



i 



(2.3) 



£ h(p) = 2J-\/l + h 2 — 2hcosp is the dispersion relation 
and K{p) is an odd function that depends on the quench 
details 



K{p) 



(ho — h) ship 



(5Jp 



+ 1 + hha — {h + hg) cosp 



Momenta are quantized as p = 2ir(n + 1/2) /L, with n 
integer. The time evolution of \*£o) is obtained by ap- 
pending the phase e l£pt to the creation operators a p . 

In RefP^l both p and p x have been computed. The 
density matrix has been identified with the GGE 

= GGE _ 
Poo — P — 

^exp^log^^np-logCl + K 2 ^))) , (2.4) 
p 



where Tip — cy.p(y.p. We remind that local integrals of mo- 
tion in Ising-like models are linear combinations of mode 
occupation numb ers^^ (analogous relations hold true 
in a quantum field theory^}, hence ensemble (2.4) is the 



GGE (1.5), exponential of a linear combination of local 
conservation laws. The diagonal ensemble p can be ob- 
tained easily by virtue of the factorization of the time 



average density matrix in the time averages of the re- 
duced density matrices (RDMs) associated to quasipar- 
ticles with opposite momenta, z.eP^l 



P^P PE 



n 

NS9p>0 



1 — n p — 
l + K 2 (p) 



bp I L—p 



(2.5) 



In RefP it was called pair ensemble for emphasizing 
the total correlation between quasiparticles with oppo- 
site momentum 

K 2 (k) 



[n k n- k ) pE = {n k ) PE = 



TTMky (^) 

We notice instead that the GGE gives a different result 
(see also RefP^l) 



{n k n_ kl 



GGE 



{n k , 



GGE 



K 2 (k) 

l + K 2 (k) 



(2.7) 



This is not unexpected since n k n^ k is a nonlocal oper- 
ator and the GGE describes only local degrees of free- 



dom (1.1). By direct comparison of Eqs ( 2.4 ) ( 2.5 1 (or 



Eqs (2.6 1 (2.7)) we see that (1.3) is not satisfied after 



quenches in Ising-like models. 

We note that in both ensembles p and p PE the von 
Neumann entropy S v .n. [p] = — Tr[plogp] is extensive. 
The entropy density of the pair ensemble is given by 



S V .N. [p 



L 



dp ( 1-K 2 (p) 
2tt \1 + K 2 (p) 



(2. 



with H{x) = 



logi±^ 



log t 



On the other 



hand the entropy density of large subsystems (of length 
I) in the generalized Gibbs ensemble (2.4) is double of 
the pair ensemble entropy density 



S v .n. [p] 



GGE] 



JO 7T VI 



K 2 {p) 



(2.9) 



£ ' J TT '"\1 + K 2 (p). 

In (the appendix of) Ref. 53 it has been pointed out that 
such discrepancies must be imputed to nonlocal degrees 
of freedom, and the authors provided an argument for the 
average of local-in-space operators to be equal in both en- 
sembles, i.e. Eq. (1.4) for quenches in the TFIC (see also 
Ref P31 for quenches to noninteracting bosonic systems). 

As emphasized in RefpS, the disagreement between 
Eq. (2.8) and Eq. (2.9) indicates that the time aver- 



age expectation values of a substantial number of ob- 
servables cannot be described through the natural finite- 
volume generalization of the GGE (which corresponds to 
remove the thermodynamic limit from Eq. (2.4) so that 



the sum runs over the momenta quantized in the finite 
volume) . It was also suggested that this could be related 
to the absence of factorization in independent uncorre- 
cted fermion-like degrees of freedom, which in the TFIC 
is attributed to the strong correlations between quasipar- 
ticles of opposite momenta (c/. Eq. ( |2.6| )). As a matter of 
fact, in the next section we show that the pair ensemble is 
factorized in noninteracting fermions, and the disagree- 
ment ( 1.3 ) is instead due to the existence of conservation 



laws independent of the quench parameter (the magnetic 
field). 
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III. PAIR ENSEMBLE 



the pair of quasiparticles wit 



As revealed by Eqs (2.51(2.6), with the time average 



1 opposite momenta emerges 



as the new elemental object. In this section we formulate 
a description in terms of "pairs" . 



First, we notice that the pair ensemble (2.5 1 describes 



nonintcracting fcrmions, indeed the eigenvalues of (2.5) 
can be written as 



v> - 



for some Vi £ [0,1], where <7j = ±1 (we notice that half 
of f's is equal to 1). However, because of Eq. (2.6), the 
nonintcracting fcrmions cannot be linear combinations 
of the fermions a that diagonalize the Hamiltonian as in 



(2.3) 



We start off with a single momentum and define two 
fermions, namely g^, which we improperly call ghost, and 
the pair fermion dp , as follows 



gl = (1 - n- p )a* p + n_ p a p 
<£ = (at + a p )aL , 



(3.1) 



with {pt , g p } = {dp, d p } = 1 and {pt, dp} = {g\, d p } = 0. 
The pair fermion represents to all intents and purposes 
the "pair" introduced qualitatively before, indeed its 
presence (absence) corresponds to the presence (absence) 
of both fermions cv± p . We have 



9p9p = n P + 
d^pdp — Ti_ p 



2n v ri- 



(3.2) 



so that the reduced density matrix at fixed momentum 
±p is given by (cf. Eq. (2.5)) 



1 



l + K 2 (p) 



+ n p n 



9p9p[ — 



(K 2 (p) - l)d%d p 



l + K 2 (p) 



which is noninteracting in gp and d p . 



In fact, the ghost and the pair fermion defined in (3.2) 



do not satisfy the correct anticommutation relations at 
different momenta: ghosts anticommute between each 
others but pair fermions commute both with ghosts and 
with other pair fermions. This problem can be easily 
solved by adding a nonlocal string to the definition of 
pair fermions: 



4= n (244-1) n ( 2 9Uk-i)di 

0<k<p 0<k^p 



(dU 



p u p 



d ] pd p ) 



(3.4) 



{9l4} = {gl,d q } = 0. 



so that 

{dp,d q } = {g ] p ,g q } = 5 P 
Finally we obtain 

Ppe = n^4exp(^logA' 2 (p)4d p -log(l+K 2 (p)) 

(3.5) 



0<k 



0<p 



which is a product of uncorrelated fermionic degrees of 
freedom. 

We also write the final Hamiltonian in terms of the 
new fcrmions 



H N s = 2e h (p)g p gl(dpd p - i 



0<p 



(3.6) 



It is rather surprising that the noninteracting represen- 
tation of p (cf. (3.5) and (2.5)) does not correspond to 



the noninteracting representation of the Hamiltonian (cf. 
Q and (lO)). 



The side effect of common conservation laws 



Comparing Eq. (3.5) with Eq. (2.4) suggests that 
Eq. (|1.3|) is not satisfied in Ising-likc models because of 



the ghost degrees of freedom. They enter Eq. (3.5 ) as the 
projector on their vacuum, and in particular 



Tr [pp E glg P ] = 0. 



(3.7) 



In addition 



9p9v 



(n p — n- p ) 2 (cf. Eq. (3.2)) is a con 



served quantity which is independent of the quench de- 



is invariant under a Bogoliubov transfor- 

with 



-n. 



(o) (o) 



= n 



n_ 



tails: n p — n_ 

mation of the a fermions, i.e. n p — n- p — ,, ., 

rip " 1 the occupation numbers of the Hamiltonian before 
the quench. The expectation value of ghost occupation 
numbers is indeed zero at any time Tr[p(t)gf.gk] = (also 
before the quench) and, despite the definition in Eq. (3.1 ) 



depends on the Hamiltonian parameters, one can even de- 
fine ghosts independent of the system details; the draw- 
back is that the Gaussian structure of pair fermions is 
lost, although the pair ensemble is still proportional to 
the projector on the ghost vacuum. In any case gj,g p is 
a conserved quantity both before and after the quench. 

does not correctly describe the expecta- 
jhost occupation numbers 



The GGE (2.4 



(3.3) ti° n value of g. 



i-'(t 



K(p) 



K 2 ( P ) 



(3.8) 



and since the right hand side is always positive, any lin- 
ear combination of gj,g p has nonzero expectation value, 
in contrast to Eq. (3.7): g p g p are intrinsically nonlocal. 



This is not totally surprising if we remind that the ini- 
tial and final Hamiltonians share half of the integrals 
of moiJOT^Hini Ghosts describe the projector on the 
eigenspaces of the common conservation laws in which 
the initial state is found; they are nonlocal because gen- 
erally the projector on an eigenspace of a local integral 
of motion is nonlocal. 

Associating zero effective temperature (A — > oo in 
Eq. ( 1.5 )) to the ghost modes identifies (3.5 1 with a gener- 



alized Gibbs ensemble. However, it cannot be produced 
as a limit of reduced density matrices, as in Eq. (1.1), 



because the limiting procedure naturally orders the in- 
tegrals of motion by their range of interaction (see also 
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RefP^) and the corresponding Lagrange multipliers do 
not have a proper thermodynamic limit. 

Similar situations can be expected also after quenches 
in interacting (integrable) models. For instance, in the 
Heisenberg XXZ chain the interaction is invariant under 
a rotation, and hence the projector of the total spin on 
the rotation axis is a local conservation law independent 
of the anisotropy parameter A of the model. Quenching 
A results in a density matrix p ( |1.2[ ) in which it is possible 
to factorize a projector. In analogy with the noninteract- 
ing case, in a basis in which the projector factorizes as 
the ghost vacuum does in Eq. (3.5 1, the remainder could 



have the form of a GGE. However, also in this case the di- 
agonal ensemble shou ld not be 'equal' to the generalized 
Gibbs ensemble ( 1.1 1. 



Another simple argument that shows that something 
can go wrong when there are common conservation laws 
is the following. If the state at late time after the quench 
is described by a GGE, the alternative time evolution 
with a local Hamiltonian that commutes with the origi- 
nal one (after the quench) and with all associated conser- 
vation laws should give rise to the same stationary state. 
However, choosing a common conservation law as new 
Hamiltonian makes the density matrix to not evolve at 
all, and hence the local properties to remain the same 
of the initial state (which, among other things, is a low- 
entangled state). 



In the next section we (re)derive Eq. (1.4) in Ising- 
like models (we formalize the argument of RefP^) and 
develop a formalism to obtain expectation values in the 
pair ensemble at C(£ _1 ). 



IV. EXPECTATION VALUES 



The pair ensemble has a very simple form when ex- 



pressed in terms of the Bogoliubov fermions (2.5) but 



computing spin correlation functions can be cumbersome 
also for the simplest operators. 

For any given observable, the time average can be 
generally replaced by the average over a finite set of 
variables, reducing the problem complexity; however 
these approaches scale exponentially with the number of 
fermions that represent the observables, i.e. with the 
distance in the two-point function of operators with a 
nonlocal fermionic representation (as the order param- 
eter in the TFIC). This makes impossible to study the 
behavior of two-point functions at large distances and 
entanglement entropies of large subsystems. 

In this section we show that working with ghosts and 
pair fermions allows us to exploit the asymptotic non- 
interacting nature of the density matrix and hence com- 
pute correlations at C(L _1 ) in polynomial computational 
time (in addition, some correlations are finally written in 
a form suitable for analytical investigations). 



A. Wick theorem in the pair ensemble: 
the removal of the ghost degrees of freedom 

We call local the observables with a spin representation 
independent of the chain length. In particular, spin cor- 
relation functions in the thermodynamic limit are local, 
regardless of the distance between operators. Strings of 
spins a? 1 ■ ■ ■ er?" are represented by strings of the Majo- 



rana fermions (1.6|, hence af and a\ are the basic local 
objects^. They are linear combinations of the Bogoli- 



(4.1) 



ubov fermions a that diagonalize the Hamiltonian (2.3) 

V p 

where 9 p is the Bogoliubov angle 

„2J(h-e*) 



£h{p) 

On the other hand, the inverse transformation of 



(3.1)(3.4) is given by 

(1 - n t)gl + n p9v 



I] (2ni-l)4e^(gl-g p ) ^ 



0<k<p 

where nt = d\d p (n 9 = gh 



n 



1), g p } 



" d 



0. 



0<fc<7T 

and in particular 



-p iL p "'p p * '"p 
The standard free fermion techniques for computing ex- 



pectation values rely on the linearity of Eq. (4.1 ), so the 



nonlinear mapping (4.2) makes calculations more diffi- 
cult. However, Eq. (4.2) is linear in ghost operators (re- 
garding e llxl ^ g as an independent operator) and, more- 
over, ghosts commute with all other pieces 



\9$ 



W n d ] 



,(t) 



. n 

Q<k<q 



1) 4 e 171 ^ 



0. (4.3) 



Thus, restricting ourselves to ghost degrees of freedom, 
we can in fact apply Wick theorem. This is the heart of 

Lemma 1 Let a, = ^ £ 0<p<7r W l (p)g p + W,(-p)gl be 

Majorana fermions (Wj(— p) — Wj(p)), linear combina- 
tions of the Bogoliubov fermions a, then 



Goa\ . . . a n Go — 



1 



L n/2 



Wx{ Pl )---W n (p n ) 



-7T<pi ,. ..,Pn <7T 

xPt([ff( R )d P( ,- R ,] < < j )Go l (4.4) 



(> 



where we indicated with [M(i,j)]j<j the skewsymmet- 
ric matrix with elements M(i,j) on the upper triangular 

part, and Go = Ilo<fc<7r 9kg\ * s ^ e P r °j ec t° r on the ghost 
vacuum. 



The proof is as follows. Since \g\, W(q)} = (cf. Eq. @ 
and Eq. (4.3)), the ghosts can be moved to the right 



preserving their order. Analogously we can move to the 
right the first projector on the ghost vacuum in Eq. (4.4 ), 



so that the ghost operators are sandwiched between two 
projectors on the ghost vacuum. We then apply the Wick 
theorem to the resulting factor, indeed 

Gog[ f) • ■ ■ g^Go = Tr g [G g{V ■ ■ ■ gW]G , 

where the trace is over the ghos t de grees of freedom, 
which results in the Pfaffian of Eq. ( 4.4 ) . Now that ghosts 



have disappeared (we are left with the projector on their 
vacuum) the operator e l7rA ^ 9 (which had just a corrective 
algebra function) can be absorbed into the definition of 
pair fermions. The removal of the ghost degrees of free- 
dom is complete. 

Lemma[T]is the connection between the noninteracting 
nature of the original model and the resulting strongly 
correlated diagonal ensemble. We are going to exploit 
Lemma [j] to prove Eq. (1.4) and then obtain the leading 
finite-size correction. The reader not interested in the 
details of the derivation can jump directly to Lemma|3]of 
Section IIVBI for the statement of the result or Section [V] 
for practical applications. 



B. Thermodynamic limit and 0(L 1 ) corrections 

If we are careful not to commute operators, we can 
depict Lemma [T] by means of Wick contractions of the 
Majorana fermions. For instance we have 



Tr^p^aiG^asa^ = Tr[pp B aia 2 a 3 a4] 



Tr[pp £ a ia' 2 a3 a' 4 ] + Tr[pp E a ia^a 3 a4 ] 



with 



Ppe = Yl 

0<p<TT 



l + {K\p)-l)a%d p 



l+K 2 (p) 



(4.5) 



the pair ensemble (3.5 ) after tracing out the ghost degrees 
of freedom. The contraction means that the sum of the 



momenta associated to the W operators ( |4.4[ ) is zero (and 
the first momentum is positive). Indeed each term of the 
Pfaffian makes the momenta equal (with opposite sign) 
two by two. We also note that each residual sum in 
Eq. (4.4) is associated to a factor 1/L. 



Definition We call tangled two contractions in which it 
is necessary to commute also the Majorana fermions of 
those very contractions in order to bring operators with 
opposite momentum close to each other (they are shown 
in Eq. (4.6)). Otherwise we say they are untangled. 



Lemma 2 The sums over the momenta of tangled con- 
tractions in Eq. (4.4) can be restricted to different mo- 



menta (in absolute value). 

The proof is as follows. In general there are two kinds of 
tangled contractions, namely 



a) 
b) 



i i 

• a, • • • a k ■ ■ ■ ai ■ 
_ i 



(4.6) 



• a, • • • a k ■ ■ ■ ar 
_ i 



which have opposite sign in Eq. (4.4), because they are 



connected by the simultaneous interchange of two dif- 
ferent rows and corresponding columns of the matrix in 
the Pfaffian (4.4), operation that changes the Pfaffian 
sign. The terms in the sums (4.4) over the momenta 



associated to the contractions (4.6) that have the same 



absolute value simplify between each other, and hence 
the sums can be constrained to run over different mo- 
menta. Since Ws with different momenta commute, we 
can reorder operators in such a way that Ws with oppo- 
site momenta are adjacent. For future reference we also 
define the tangled group as follows: 

Definition We call tangled group of a contraction the 
set of contractions in which momenta are forced to be 
different (by Lemma [2]) from the momentum of the con- 
traction. 

In terms of ghosts and pair fermions the Majorana 
fermions (1.6) can be written as (cf. Eq. (4.1)) 



-i= £ Afy\p) gp + Afy\-p)gl. (4.7) 

VL 0< P <7T 



with (p > 0) 



Af(p) = co S (p£ - f )[1 - iif] + i sm(p£ - -f)[ T ; + t; 
A»(p) = sm(p£ + |)[1 + ir»] -icos(p£ + \)[ 



2 / \T%> \ 



(4.8) 

and Af v) (-p) = (Af v) )Hp). The Pauli matrices t*' v < z 
are given by ([t*,t|] = 2i5 pq ei jk T^) 



0<fe<p 



t c inN„ 



r z p = 2n a p 



1 . 



(4.9) 

We agree to order the Majorana fermions in the ob- 
servables a? 1 • • • a?" as follows: 



a j < a )+i < a r+i' < 'V 



(4.10) 



where j,j',£,£' > 0. In this way some simplifications, 
which are otherwise difficult to recognize, are evident. 
The contraction of adjacent Majorana fermions reads as 



7 



with 



0<p 



L 



0<p 



(4.12) 



where the apostrophe in the sums is used to remind that 
the momentum is different from any other momentum in 
the tangled group. We also write the contraction (which 
now means equality of momenta) between / and g oper- 
ators: 

2 - s „ 2 „ _ 2 - 

fefi' = -fft+l' j Jlffi' = y9i'-i ; ffiffi' = TJ£'-t > 

1-1 " 1-1 " 1-1 (4.13) 

where the tangled groups of the resulting / and g opera- 
tors are the union of the tangled groups of the contracted 
operators. 



We now prove Eq. ( 1.4 1. The operator norm of / and g 



contraction gives a contribution 0(L x ) (cf. Eq. (4.13 1). 



Spin correlation functions are expectation values of the 
product of a finite number of Majorana fcrmions, hence 
the number of contractions that come out of Eq. (4.4) is 
independent of L. At the leading order O(L ) we can 
neglect such contributions and assume that all Majorana 
fermions have distinct momenta (we say that they have 
"maximal tangled group" ) . The average of products of / 
and g factorizes in the product of the averages, in which 
sums run over distinct momenta: 



0<p 

2 W ,„ . „ , 1 -K 2 (p) 



Tr[p PE ge} =-£ ^ caa(£p + O p )- 



(4.14) 



0<p 



however in the thermodynamic limit the constraints dis- 
appear 



" dp 



o+iQp 



l - K 2 (p) 
TTkHp) 



(4.15) 



since the number of sums and terms of the Wick ex- 



pansion (4.4) is finite and independent of L. Eq. (4.4) 



and Eq. (4.15) are nothing but the Wick decomposition 
in the GGE (12. 41, that is to say Eq. il.ty for noninter- 



acting Ising-like models; but this is not the end of the 
story. Eq. (4.14) is in fact sufficient to characterize the 
C(L _1 ) correction too (see Appendix |A| and, in addition, 
the terms contributing at 0{L^ r ) have no more than two 

/'s. 



We now rewrite Eq. (4.14) in a more practical way. 



The terms of the Wick expansion without / operators are 



trivially described by the correlators (4.14), also without 



the constraints on the momenta, because with the order- 
ing (4.10) g operators are tangled (see Appendix |A|) . On 



the other hand, the terms with two /'s have the following 



form 



(ftfe 



4 

L 2 



E 

0<fe,fe / 

k^k' 



cos(£fc) cos(£'fc') 



/ pi^k,k' 



where we had to exclude the momenta associated to the 
/'s in the remaining expectation value (of g operators). 
At the first order, the two momenta in the expectation 
value are decoupled 

(• • 'Wm' = (■••) + (QW + Q(k'))/L + 0(L- 2 ) 

where Q is some bounded function, which could be writ- 
ten explicitly by series expanding the Pfaffian in (■••). 
In fact we don't need the exact expression of Q, because 
the leading correction to (■ • •) is C(L _1 ), but the sum 
over the other momentum is 0{L^ r ) as well, resulting in 
(D(L~ 2 ) contribution. Finally we obtain (£,£' > 0) 



(4.12) is less then or equal to 1. On the other hand, each (fifi* ■ ■ ■) 



; cos(£fc) COs(f fc) (• • •) 



26, 



L 2 ^ 

0<k 

(4.16) 

As a matter of fact, it is convenient to represent the Kro- 
necker delta as a sum, like in the central expression of 
Eq. (4.16). However, there is no need to sum over L/2 



elements, indeed we can write 

(ftft> •■■>~ -yjH cos ( fk ) cos ( rt ) <•••>, (4-17) 

0<k 

for any L > I + £' and the sum is over L/2 momenta 
k = 2ix(n + 1/2)/ L, with n integer. Considering local 
observables, the contractions of a x and a y fermions give 
rise to operators fe in which the maximal £, let us say 
£m, is independent of L; it is sufficient to choose L > 2£m 



to have a representation (4.17) valid for all terms of the 



Wick expansion (4.4) and for any chain length L. In this 



way we can easily isolate a factor 1/L from the expres- 
sions. 

We note that each term of the (averaged) sum in 
Eq. (4.17) is a term of the Wick expansion of a non- 



interacting model with modified correlators, and taking 
into account the algebra of operators (Eq. (4.17) is valid 
only for operators ordered as in Eq. (4.10)) we finally 
obtain: 

Lemma 3 At C(i _1 ) the expectation value of local ob- 
servables in the pair ensemble can be written as the aver- 
age of the expectation values computed in L/2 Gaussian 
states: 



{ai ■ ■ ■ a 



n/PE 



7 ^(a 1 ---a n )^+ (i- 1 ), (4.18) 



0<k 



where (sgn(0) = 0) 



(a*+X)j = ( aV r a r+i) k = &o + i\l ^sgn(£) cos{k£) 



* dp itM l-K\p) 
^2tt l + K 2 (p) 

(4.19) 



and L is larger than the maximal difference between the 
indices of the a x ^ fermions in Eq. (4.18). 



The leading finite-size correction ©(L" 1 ) can be ob- 
tained by series expanding at the second order the Pfaf- 
fian associated to the expectation values at fixed mo- 
mentum k, as we are going to do for the longitudinal 
two-point function and the Renyi entropy. 



V. FINITE-SIZE CORRECTIONS 

We consider first some operators with a local fcrmionic 
representation. 

Observables that are quadratic in the Majorana 
fermions ( 1.6 ) 



v v 



etc. 



have exponentially small finite-size corrections, which 
come out of the substitution of sums with integrals in 
Eq. ( 4.14[ ), being e lSp and K 2 {p) smooth periodic func- 
tions (at least for quenches between noncritical models). 
These are the same corrections that arise considering the 



finite volume generalization of Eq. (2.4 1, so they are al- 
ways present but independent of the time average. // 
there were no common conservation laws, we would ex- 
pect only this type of corrections] therefore, we under- 
stand any additional correction as the effect of conserva- 
tion laws independent of the quench parameter. 

We now consider operators consisting of four Majo- 
rana fermions. In this simple case we can compute the 
expectation value explicitly. We obtain 



(a v ja v J+e a*, +e 



' PE 



GGE 



(5.1) 



Ignoring the exponentially small correction discussed 
above the finite-size correction is nonzero only if I = £', 
i.e. for the two-point function of the operators o~ x o~ x , 



and o\a\. 



A. Transverse correlation 

The two-point function of the transverse field is the 
expectation value of four fermions: 



hence from Eq. (5.1) we have 



CJ>eW ~ Cq GE {C) ~ 2/L . 



(5.2) 



(5.3) 



The difference between the 2-point transverse correla- 
tions in the GGE and in the pair ensemble is independent 
of the distance. It is worth to compare this simple result 
with the approach to the stationary state. We remind 



that in the GGE the connected two-point function of the 
transverse field decays exponentially with the distance^. 
In RefPS it has been observed that in order to extract 
the correlation length £ z of the transverse field one has 
to wait a time exponentially larger than the distance. In 
finite systems, to avoid quantum revivals, the time must 
be smaller than the chain length L and hence log L ^> I. 
With the time average we obtain an analogous behavior. 
Indeed, because of Eq. (5.3), the exponential decay can 
be 'seen' only if e~ 



L ' 



that is to say 



< log L . 



(5.4) 



B. Order parameter 2-point function 

In this section we consider the longitudinal two point 
function 



C x (£) 



(5.5) 



In contrast to the transverse correlation the leading 
finite-size correction does depend on the distance, e.g. 



'l u 2l pe 
'1 a 3/PE 



'\ U 1IGGE 

x x\ 
T l"3/ GGE 



2 
L 



(5.6) 



Computing C x {tj is much more complicated than com- 
puting C z (£), since the longitudinal correlation is non- 
local in terms of the Jordan- Wigner fermions. We use 
Lemma [3] to compute the leading correction. In each 
Gaussian state of eqs (4.18 )(4.19b the two point function 



has the pfaffian representation^ 



/ * X * X \ 



Pf 



with 

Gmn Tro— n— 1 — 



71 dp 



, G t F (k) 



i(m—n — l)p-\-i9 p 



(5.7) 



l-K\p) 
1 + K 2 (p) 



i—n — * 
m, n = 1 



sgn(m — n) cos(fc(m — ri)) 



J. 



The function 

X 

characterizes the leading correction to Cp E 
C x PE {t) = C&geW (l + ^ + 0(L-i 



= - lim L 2 ^logC^) 



(5.8) 
(5.9) 

(5.10) 



It can be computed by series expanding the Pfaffians in 
Eq. (4.18) at the second order in L -1 / 2 and averaging 
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1x10 



1x10' 



r 



1x10 




FIG. 1. The coefficient of the leading finite-size correc- 
tion for the time average longitudinal two-point correlation 
for various quenches (ho,h) in the TFIC. In quenches from 
the ordered phase xf displays a clean power law behavior, 
whereas in quenches within the disordered phase (the black 
circles represent — xf) it grows exponentially with the dis- 
tance. Quenches across the critical point from the disordered 
phase (green down-pointing triangles) are characterized by os- 
cillations, however we recognize a parabolic growth for almost 
all lengths. 



C. Renyi entropy Sz 

In this section we compute the leading finite-size cor- 
rection for the Renyi entropy S 2 = — \ogTr[p 2 }, where pi 
is the RDM of t adjacent spins. We use that at the lowest 
order in L^ 1 the RDM can be written as (c/. Lemma [3]) 



1j 

0<k 



(5.13) 



where p[r( fc )] is the Gaussian matrix with correla- 
tions ( |4.19[ ). In particular, is the block Toeplitz 



matrix 



r (fe) = 
ij 



Pi-j 

-Hj-i ~<f>i-j 



i,j = !,...,£, (5.14) 



with 7 and <j> of Eq. ( 5.8 ) . We need to compute the second 
moment of the reduced density matrix Ti[(pf E ) 2 ] 



Tr[(pn 21 



L 2 



£ Tr[p[rW] P [r^')]] . 



(5.15) 



0<k,k' 



The trace of the product of two Gaussian density matri- 
ces is given bj^21 



over the auxiliary momentum. We finally obtain 



0<k 



E ,k+j ~ [@ X ]kl ■ (5-11) 

i,j,k,l= 1 



In Fig. [T] we report the numerical data for quenches 
within and between the ferromagnetic and paramagnetic 
phases of the TFIC. The leading correction can be ne- 
glected in the following regimes: 



'l 1 ' 2 h<l 

: < i 1 / 3 h <l,h>l 
logL ho,h > 1 . 



(5.12) 



Tr[p[r]p[f]] = det 



1 + rr 



1/2 



(5.16) 



so the leading finite-size correction can be computed by 
series expanding the determinants in Eq. (5.15) treating 
(/>'s as perturbations. The function 



02U) = lim L 



dL 



PE 



characterizes the leading correction to S 2 
S? E W = S§ GE {i) - ^ + 0{L- 2 ) 



After some algebra we obtain 

1 



0<fc 



R 



-1 



R f J R- 1 -+ 



R 



(5.17) 



(5.18) 



-F ( ~ k) ] = 
t ' 



In Ref d 51 * 53 ^ it has been shown that after quenches origi- 
nating in the ferromagnetic phase the stationary value of 
C x (£) emerges after a time that scales as a power law 
with the distance £; on the other hand the timescale 
after which the stationary behavior reveals itself after 
quenches within the paramagnetic phase is exponentially 
large. The finite size corrections (5.12| display the same 



behavior, providing further indications that large finite- 
size corrections in finite chains and slow relaxation in the 
thermodynamic limit could have the same root, i.e. in this 
specific case the existence of common conservation laws 
before and after the quench. 



0>i+t,k+o - ,k+i 

i,j,k,l 



R- 



R* 



L R- 



R 



ki 
(5.19) 



with F the matrix defined in Eq. ( 5.8 ) and Rij — ji-j . In 
Fig. [2] we report the numerical data for quenches within 
and between the ferromagnetic and paramagnetic phases. 
The leading correction can be neglected when I <C L 1//2 , 
but in comparison with the extensive part of the Renyi 
entropy we have 



<L. 



(5.20) 



10 
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FIG. 2. The coefficient of the leading finite-size correction for 
the time average Renyi entropy S2 for various quenches (ho, h) 
in the TFIC. In all cases the leading correction oi(£) ~ t 2 . 



the system size as a multiplicative factor (which in Ising- 
like models is 1/L) and then take the limit of infinite 
subsystem in the infinite system. In this limit p GGE is 
substituted with p , which is the exponential of the 
local conservation laws. However, by construction, the 
two ensembles share the same local properties, so they 
must give the same expectation values of local integrals 
of motion Ij, i.e. (cf. Eq. (1.5 1) 



AS vN = Trpog^ Sp] + 0(5pj) 

Tv[\ogp GGE Sp] = Tr[Q^ \, J I j )(p GGE 
3 



p)] = 0. 



(5.24) 

Thus we expect that either the first correction drops off 
to zero for large lengths or it is exactly zero at any length, 
as in the TFIC. However this argument can be used only 
to exclude the leading power law correction in L and, 
e.g., does not apply when all (local) observables have 
exponentially small finite-size corrections. 



Eq. (5.20) is the time average analogous of the large time 
behavior of the Renyi entropy^ S 2 (£,t) ~ S 2 {£,oo) + 
0(£ 4 /t 3 ), where the correction can be neglected only if 
I <C Jt. We point out that the leading finite-size correc- 
tion is negative, as one could expect having in mind the 
entropy of the whole chain (cf. eqs ( 2.8 ) ( 2.9 1 for the von 
Neumann entropy). 



D. Von Neumann entropy 

In this section we show that the correction for the von 
Neumann entropy is (D(L~ 2 ). We indicate with Spe the 
difference between the RDMs in the pair ensemble and 
in the GGE 



APE) 



(GGE) , 

Pe + ope 



(5.21) 



At the first order in Spe (which corresponds to the first 
order in 1/L at fixed £ and large L) we have 

Svn\p¥ E) ] = S vN [pf GE) \ - Tr[logpf GE 6 Pe ] + 0(L- 2 ) . 

(5.22) 

The RDM in the GGE is the exponential of a two-form, 



so the trace in Eq. (5.22) is nonzero only if Spe has some 
quadratic contribution. This is however impossible, be- 
cause the expectation value of fermionic quadratic oper- 
ators is the same in the two ensembles. Thus we obtain 



{£) 



qGGE 



(£)+0(L- 2 ). 



(5.23) 



There are indications that this result can be partially 
generalized to generic quenches in which Eq. (1.4) is sat- 



isfied and the local properties of the system at late times 
after the quench are described by a GGE. 

The argument is the following. At the first order we 
can extract from Spe = p 



VI. CONCLUSIONS 

In this work we have shown that quenches of global 
parameters in periodic noninteracting chains are gener- 
ally characterized by (infinite) conservation laws com- 
mon to the Hamiltonians before and after the quench. 
We provided evidence that they have an impact on the 
finite-size corrections for time average expectation val- 
ues. Considering in particular quenches in the quantum 
Ising model, wc identified the regimes in which the finite- 
size corrections are negligible and found that large cor- 
rections correspond to slow relaxation dynamics in the 
infinite chairPM This indicates that in some cases slow 
relaxation could be related to the fact that the initial 
state is (close to be) eigenstate of local conservation laws 
of the final Hamiltonian. 

In noninteracting models the slowing down could be 
amplified by the fact that there are infinite common con- 
servation laws, so this issue need further investigations. 
To this aim, e.g., it could be worth analyzing quenches 
in which the noninteracting representation of the Hamil- 
tonian depends on the quench parameter (in contrast to 
the Jordan- Wigner transformation (1.6)). 



Finally, we did not calculated either the time average 
two-point functions or Renyi entropies when the typical 
length is order the chain length. Computing the (full) 
finite size scaling appears to be a very hard problem. In 
fact, semiclassical theories turned out to reproduce the 
asymptotic behavior of entanglement ent ropies a nd corre- 
lation functions after a quantum quencfP- 1 * 47 * 68 * ^; since 
there are sizable finite-size corrections for the time avcr- 



pf GE the dependence on 



age entanglement entropy density (cf. Eqs (2.8)(2.9|) wc 
believe that a semiclassical approach might be successful 
also in this context. 
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Appendix A 



Here we prove that the factorization (4.14) holds true 
at 0(L^ r ). First we observe that in each term of the 



Wick expansion of the operators ordered as in (4.10) all 
g's are tangled, indeed untangled <?'s appear only if a y 
fermions come both before and after a x fermions, condi- 
tion never satisfied if Majorana fermions are ordered as 
in (4.10). Thus, any contraction involves / operators. In 



addition, from Eq. (4.13) we see that the leading correc- 
tion involves no more than one contraction. We 
also note that the observables with nonzero expectation 
values have an even number of f. This is a trivial conse- 
quence of the fact that, whatever contractions are done, 
the contribution from / and g operators ( 4.12[ ) is always 
real. The adjoint of a string of 2n distinct Majorana 
fermions is the same string multiplied by (— l) 1 ^ 2 ™" 1 ); 
on the other hand each g operator is associated with a i 
(cf. Eq. (4.11)), and hence the expectation value takes 
a factor (— l) n ~ n f , where n/ is the number of / opera- 
tors. In order to have a nonzero expectation value the 
two factors have to be equal, i.e. nf must be even, as 
required. 

Now we show that in the terms contributing at C(L _1 ) 
no more than two f 's are present. 

First we estimate the contribution from /'s with maxi- 
mal tangled group. This can be inferred from the relation 



{t} n ee (I)" E cos (^) ^ 0(L-L(«+D/ 2 J) , (Al) 

0<P1 , • - • ,P2n 
Pi^Vj > k l >■■■ k m 

valid for any n < L, m < L, ^ > and momenta hi. 



Eq. ( Al ) can be proved writing the recurrence equation 



{£} n = 5^(1 + |(1 - n- m)) {£}„_! - 

m n— 1 

[2]>>s(^)mn-i + E E m-^' 



1 - <^„0 



L 



X 



j = l cr = ±l 



and noticing that we don't get a factor j only when the 
Kronecker delta is 1; since £j ^ the delta can be 1 only 
for the terms like the last one, where lengths are shifted. 
In the worst case we can sum over [(n + 1)/2J momenta 
without satisfying the delta, ending up with a term in 
which all remaining deltas are equal to 1. Thus, we have 



Eq. ( Al ) and in particular we see that both a single / 
(n = 1) and two /'s (n = 2) with maximal tangled group 
contribute at ©(L" 1 ), whereas more than two /'s give 
subsleading contributions. 

On the other hand, if /'s have not maximal tan- 
gled group they are contracted. Each contraction gives 
0(L~ r ) contribution (at best, cf. Eq. (4.13)), hence we 
conclude that, in any case, more then two /'s correspond 
to subleading terms. 

Since g's are tangled, terms without /'s are simply de- 
scribed by Eq. (4.14). We focus on terms with two j's. 



The case of one / contracted with a g and one / with 
a maximal tangled group is subleading, since both oper- 
ations are 0(L~ r ). The contraction of the / operators 
is subleading as well because the ordering (4.10) makes 



positive the index of /'s, forbidding the appearance of fo 
as a result of the contraction (cf. Eq. (4.13)), which is 
the only term that would give contribution. 

In conclusion, the tangled group of both f's must be 
maximal, so Eq. (4.14) is in fact sufficient to characterize 
the 0(L _1 ) correction too. 
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